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ABSTRACT 



Context. The use of SN rates to probe explosion scenarios and to trace the cosmic star formation history received a boost from a 
number of synoptic surveys. There has been a recent claim of a mismatch by a factor of two between star formation and core collapse 
SN rates, and different explanations have been proposed for this discrepancy. 

Aims. We attempted an independent test of the relation between star formation and supernova rates in the extreme environment of 
starburst galaxies, where both star formation and extinction are extremely high. 

Methods. To this aim we conducted an infrared supernova search in a sample of local starburts galaxies. The rational to search in 
the infrared is to reduce the bias due to extinction, which is one of the putative reasons for the observed discrepancy between star 
formation and supernova rates. To evaluate the outcome of the search we developed a MonteCarlo simulation tool that is used to 
predict the number and properties of the expected supernovae based on the search characteristics and the current understanding of 
starburst galaxies and supernovae. 

Results. During the search we discovered 6 supernovae (4 with spectroscopic classification) which is in excellent agreement with the 
prediction of the MonteCarlo simulation tool that is, on average, 5.3 ± 2.3 events. 

Conclusions. The number of supernovae detected in starburst galaxies is consistent with that predicted from their high star formation 
rate when we recognize that a major fraction (~ 60%) of the events remains hidden in the unaccessible, high density nuclear regions 
because of a combination of reduced search efficiency and high extinction. 

Key words. Stars: supernovae: general - Galaxies: starburst - Galaxies: star formation - Infrared: galaxies - Infrared: stars 



1. Introduction 

The rate of supernovae (SNe) is a key quantity in astrophysics 
that provides a crucial test for stellar evolution theory and an 
input for the modeling of galaxy evolution with direct impact 
on the chemical enrichment and the feedback mechanism. Core- 
collapse SNe (SN CC), because of their short-lived progenitors, 
trace the current star formation rate (SFR). Conversely, for an 
adopted SFR, measurements of the SN CC rates give informa- 
tion on the mass range of their progenitors as well as the slope 
of the initial mass function at the high mass end. SN la, resulting 
from the thermonuclear explosion of a white dwarf in a binary 
system, show a wide range of delay times from star formation to 
explosion. Therefore, the SN la rate reflects the long-term star 
formation history of the parent stellar system. Recently, it has 
been claimed that a significant fraction o f SN la have a short de - 
lay time, possibly as short as 10 7 years dMannucci et al.l l2006). 
Like for CC SN, the rate of such prompt SN la events is expected 
to be proportional to the current SFR. 

In one of the early a ttempts to compare the SN and SF rates, 
ICappellaro et alJ (Il999l) found that the SN CC rate in galaxies 
with different U — V color matches the predicted SFR when 



* ESO proposal: 083.D-0259, 085.D-0335, 085.D-0348, 087.D- 
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adopting a mass range 10 M Q < M < 40 M for the SN CC 
progenitors. 

In the last decade there was a enormous improvement in the 
measurement of the cosmic SFR with the careful co mbination 
of many different probes (eg. Hop kins & Bea com 2006). A most 
relevant feature is that the SFR reaches a maximum at a redshift 
Z ~ 1 and hereafter begins to decrease down to the current rate 
which is over one order of magnitude lower than at peak. 

A significant effort was also devoted to the measurement 
of the cosmic SN rate: although much of the focus was 
for type SN la, a few estimates of the SN CC rate s were 
also published both for the local Universe dLi et al. 
and at high redshifts (|Dahlen et all 12004]: ICa ppellaro 
[ 20051; IBotticella et alJ2008tlBazin et alj2009l: lGra ur et al.ll201 
llvlelinder et al.lr2012tlDahlen et al.ll2012l) . While the new local 
SN CC rate confirms previous results, with a much better statis- 
tics and lower systematic errors, the evolution with redshift was 
found to track very well the SFR evolution, considering the large 
uncertainties in the extinction corrections. Again, to best match 
the observed SN and SF rates it was argu ed that the lower limit 
for SN CC progenitor h ad to be ~ 10 M Q ( IBotticella et al.ll2008t 
iBlanc & Greggiol2008l) . 

At about the same time, following a different line of research, 
the analysis of archival images allowed the identification of the 
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precursors for a number of nearby SN CC. From the often very 
scanty but precious photometry, and using stellar evolution mod- 
els, one can estimate the SN precursor mass. The uncertain- 
ties are in general quite large, as confirmed from the discrep- 
ancy in the mass estimates from different groups, but this anal- 
ysis suggests a lower limit for SN CC progenitors of 8 + 1 M 
(Smartt 2009). If this value is adopted, the observed SN rates 
would result a factor two smaller than those expected from the 
observed SFR. Thi s was identified by s ome authors as a "SN 
rate problem" (e.g. iHoriuchi et al.ll201ll While one should re- 
mind that the uncertainties on SFR rate calibratio ns are still large 
(Bott icella et al.l2012trKennicutt & Evansl2012l) . it also true that 
there is a number of possible biases in the SN rate estimates. The 
two most severe are the possible underestimate of a large popu- 
lation of faint SN CC and/or the underestimate of the c orrection 
for extinction (IHoriuchi et al.ll20lTt iMattila et al.ll2012l) 

In particular, IMannucci et al l (120071) : ICresci et ail (120071) 
and, more recently, Mattila et al. (2012) argued that a significant 
fraction of SN CC remains hidden in the nuclear region of star- 
burst galaxies, with a loss of up to ~70-90% in the highly dust- 
enshrouded environments of (ultra-)luminous infrared gal ax - 
ies( U/LIRGs). This effect is expected to be more important at 
high redshift because of the larger fraction of starburst galaxies. 
Indeed, when a correction for this hidden SN fraction is included 
in the rate calculation the discrepancy between SN and SF 
rates at high redshi fts seems to disappear (Meli nder et al.ll2012t 
iDahlen et alJl2012t the "missing frac t ion" c orrection adopted in 
these works was from IMattila et al.l (120121) . It is currently un- 
clear if this effect is large enough to explain also the discrepancy 
observed in the local Universe with somewhat conflicting evi- 
dences from the statistics of SNe in the Local Group galaxies 
dBotticella et alj|2012t IMattila et ai1l2012l) and large sample SN 
searches ( Li et alJ201 la ). 

Entering in this debate, we planned for an infrared SN search 
in a sample of local starburst galaxies (SBs). The idea was to 
verify the link between SN and SF rates in an environment where 
star formation is very high, 1-2 order of magnitude higher than 
in normal star-forming galaxies. By observing in the K-band we 
were aiming to reduce the bias due to extinction (Ak ~ 0.1 Ay). 

The idea is not new. A first attempt of a dedi cated SN search 
in SB s was performed in the optical band by Rich mond et al.l 
(1998). During the search only a handful of events were de- 
tected leading the authors to conclude that the rate of (unob- 
scured) SNe in SBs is the sam e as in quiescent galaxies. A sim- 
ilar conclusion was reached bv lNavasardvan et al.l (12001 1) , again 
based on optical dat a. As for infrared SN search, after a few un - 
successful attempts dGrossan et al.ll 1999t iBregman et al.l l2000). 
the first results of a sy ste matic search in SBs wer e reported by 
iMaiolino et all (120021) and IMannucci et alJ (120031) . They found 
that the observed SN rate in SBs was indeed one order of mag- 
nitude higher then expected for the galaxy blue luminosities but 
still 3-10 times lower than would be expected from the far in- 
frared (FIR) luminosity. Among the possible explanation for the 
remaining discrepancy, they suggested extreme extinction in the 
galaxy nuclear regions (Ay > 25mag), which would dim SNe 
even in the near-IR, and insufficient spatial resolution to probe 
the very nuclear regions. The reliability of the use of NIR search 
for obscured SNe in the nuclear and circumnu clear regions of ac- 
tive st arburst galaxies was also investigated bv lMattila & Meikla 
J2001 taking into account in particular the problem of extinc- 
tion. They conclude that with a modest investment of observa- 
tional time it may be possible to discover a number of nuclear 
SNe. A negative search for transients in NICMOS images re- 
trieved from the Hubble Space Telescope archive suggests that 



the same biases like ly affect also space -based, high spatial reso- 
lution observations ( Cresci et al.112007 ). 

The same approach was used by Mat tila et al.l {2007b) but 
with ground based, adaptive optics (AO) assisted observations. 
The applic ation of this technique led t o the discovery of a hand- 
ful of SNe dKankare et alj|2008ll2012l) but not yet to an estimate 
oftheSNCCrate. 

Until now, about a dozen SNe have been discovered by IR 
SN searches, not all with spectroscopic confirmation. The num- 
ber is higher if we include also events first detected in the optical 
and re-discovered by the IR searches. Therefore the statistics is 
still very low and many of the original questions are still unan- 
swered. This gave us the motivations to make a new attempt ex- 
ploiting the opportunity offered by HAWK-I, the infrared camera 
mounted at the ESO VLT telescope. 

The paper is divided in two parts: the first part describe the 
observing program, namely the galaxy sample and the search 
strategy in Sect. 12.11 the data reduction in Sect. 12.31 the SN 
discoveries and classification in Sect. 12.41 while in Sect. 12.51 we 
detail the procedure to estimate the search detection efficiency. 
The second part is devoted to the description of a simulation 
tool which is used to predict, based on our current knowledge of 
SBs properties and on the specific features of our SN search, the 
number of expected SN detections (Sect. [3]!. Finally, we compare 
the number and properties of the expected and observed events 
(Sect.|4]i and draw our conclusions (Sect. [5}. 

Throughout this paper we assume the following cosmolog- 
ical parameters: Hq = 72kms~ 1 Mpc -1 , Oa = 0.73 and Qm = 
0.27. 



2. The SN search program 

2.1. Galaxy sample 

Starbursts are galaxies with very high star formation rate, of the 
order of 10-100 Mq yr~' compared to the few Mq yr~' of nor- 
mal star forming galaxies in the local universe. Given that in a 
typical galaxy the very high SFR will rapidly consume the gas 
reservoir, it is thought that the starburst is a temporary phase in 
the galaxy evolution. The fact that many SBs are in close pairs or 
have disturbed morphologies point to the interaction as a dom- 
inant, alth ough possibly not unique, reason of the phenomena 
(Gallagherl ll993l) . The ultra-violet radiation from young, mas- 
sive stars heats the surrounding dust and is re-emitted in the far 
infrared. Indeed the most luminous SBs in the local Universe 
are LIRGs , with 1 1 < log(L/j;/Lm) < 12, and ULIRGs, with 
log(L /s /L ) > 12 (ISanders & Mirabelll 19961) . 

For our projec t we selected from the IRAS Revised Bright 
Galaxy Sample (ISanders et al.ll2Q03h a sample of SBs with to- 
tal infrared (TIR) luminosity log(L7-/«/Lo) > 11 and redshift 
z < 0.07. With the additional requirement that the targets are ac- 
cessible from Paranal in the April to September observing sea- 
son (to fit in one of the ESO allocation period) we retrieved a 
sample of 30 SBs. 

The list of SBs is reported in Tab. Q] Along with the galaxy 
name and equatorial coordinates (cols. 1-3) we report the helio- 
centric redshift (col. 4), log Ljir and log Lb (cols. 5 and 6; cf. 
Sect. 13.1. Il l, the Hubble type (col. 7), the SFR and the expected 
SN rates (cols 8, 9) derived from Ljm as described in Sect. l3.1.T1 
Galaxy data have been retrieved from NETjfl In the last column 



1 The NASA/IPAC Extragalactic Database (NED) is operated by the 
Jet Propulsion Laboratory, California Institute of Technology, under 
contract with the National Aeronautics and Space Administration. 
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Fig. 1: Distribution of the B and FIR luminosities for the SB 
galaxies of our sample. 



we listed (in boldface) the designation of the SNe discovered in 
our search which are the basis for our analysis. For completeness 
we also list (in italics) the SNe discovered by other SN searches 
outside our monitoring period. The distribution of Lb and Ljm 
are compared in Fig. Q] showing that, as typical for SBs, Ljm is 
on average a factor ten higher than Lb, whereas for normal star 
forming galaxies Ljm ~ Lb- We notice that almost all galax- 
ies are LIRGs and only two are ULIRGS. Most galaxies of the 
sample are isolated (~ 60 - 70%) while the remaining are dou- 
ble/interacting galaxies or contain double nuclei, signature of a 
recent merger. Several galaxies of the sample are asymmetrical, 
disturbed, or show warps, bars and tidal tails. 

2.2. Search strategy 

To search SNe in the selected SB sample we used the HAWK-I 
instrument installed at the ESO VLT telescope at Cerro Paranal 
(Chile). HAWK-I is a NIR (0.85 - 2.5 /mi ) wide-field imager 
with a mosaic of four Hawaii-2RG detectors. The total field of 
view is 7.5' x 7.5' with a scale of 0.106 "/pix. Even in poor 
seeing conditions (> 1.5 arcsec) the instrument allows to achieve 
S/N ~ 10 for a K=20 magnitude star with a 15 min exposure. 

The infrared light curves of SNe evolve relatively slowly, 
remaining within one/two magnitudes from maximum for 
two/three months (Mattila & Meiklel l200ll) and therefore an IR 
SN search does not require frequent monitoring. We planned for 
an average of three visits per galaxy per semester, for a total of 
80-100 visits. The monitoring campaign was scheduled in ser- 
vice mode and we did not set tight constraints for the sky con- 
ditions. This and the relatively short duration of the observing 
blocks made the program well suited as filler. We notice that we 
had no influence on the actual scheduling of the observations 
which followed the rules of the ESO service mode scheduler. 

Eventually, the fraction of useful observing time was 100% 
of the allocated time in the first season, and 70% in the second 



and third semesters. The log of the observations is reported in 
Tab. |2] where for each galaxy we list the epoch of observations 
(MJD), the seeing (FWHM in arcsec), and the minimum and 
maximum magnitude limit for SN detection across the image 
(cf. Sec. 12.5b . In total, we obtained 210 K-band exposures (ex- 
posure time 15min), with an average of about 3 visits per galaxy 
per semester. Because of the time loss, three galaxies were not 
monitored in the last two seasons. 

It turned out that the average image quality was quite good: 
for ~ 90 % of the exposures the seeing was less than 1.0", with 
an average FWHM across the whole program of 0.6". 

2.3. Data reduction and analysis 

For data reduction and mining of the HAWK-I mosaic images we 
developed a custom pipeline that integrates different, publicly 
available, recipes and tools in a Python environment. 
The pipeline consists of four sections: 

1. pre-reduction, astrometric calibration and production of 
the stacked mosaic image. For these steps we use the 
ESO HAWK-I pipeline recipes in EsoRex, the ESO Recipe 
Execution Tooq 

2. subtraction o f images taken at different epochs using ISIS 
(lAlardll20 00) for the PSF matching; 

3. search for t ransient candidates in the difference image using 
Sextractor dBertin & Arnoutsl Il996l) . The candidates were 
ranked based on their Sextractor measured parameters and 
submitted to the operator for visual inspection and valida- 
tion; 

4. estimate of the detection efficiency through artificial star ex- 
periments performed for each of the search images (details 
in SectlOTl. 

The raw images were retrieved from the ESO archive as soon 
as they became available, and immediately reduced to allow for 
activation of follow-up spectroscopy of transient candidates. 

For the pre-reduction, we followed the reduction cascade de- 
scribed in the HAWK-I pipeline manual including dark subtrac- 
tion, flat field and illumination corrections, background subtrac- 
tion, distortion correction, astrometric offset refinement, com- 
bination of the different exposures and stitch of the 4 detec- 
tors in a single mosaic image. Actually, it turned out that the 
ESO pipeline recipes for background subtraction and offset re- 
finement do not provide satisfactory results for our images. The 
main reason is the extended size of our sources and the conse- 
quent large dithering we had adopted. To address this issue we 
implemented custom recipes for the two afore mentioned reduc- 
tion steps. 

The most critical step of the data reduction is the image sub- 
traction, in particular in the proximity of the nuclear regions 
of the galaxies. First of all we need to choose a proper refer- 
ence image, usually the image with the best seeing obtained at 
least three month before (or in some case after) the image to be 
searched. We also need to choo se the proper paramet ers for the 
image difference procedure (see iMelinder et aO 2012 for an ex- 
tensive discussion). An additional problems arises because in the 
distributed version of ISIS, the program automatically selects the 
reference sources for the computation of the convolution kernel. 
Owing to the small number of sources in our extragalactic fields, 

2 http://www.eso.org/sci/software/cpl/esorex.html 

3 ftp://ftp.eso.org/pub/dfs/pipelines/hawki/hawki-pipeline-manual- 
L8.pdf 
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Table 1: The SB galaxy sample. The last column report our 6 SNe (in bold) with other SNe discovered previously in the galaxy 
sample. 



galaxy 
designation 


R.A. Dec. 
J2000.0 


redshift 


log L TIR 

Pol 


log Lb 
[L©] 


Hubble 

type 


SFR 

[M yr'] 


SN rate 
[SNyr 1 ] 


SNe 


CGCG01 1-076 


11 21 13.3 


-02 59 08 


0.025 


11.28 


10.35 


2.9 


32.1 


0.38 




CGCG043-099 


13 01 49.9 


+04 20 01 


0.037 


11.59 


10.51 


3.4 


65.4 


0.77 




ESO148-IG002 


23 15 46.6 


-59 03 14 


0.045 


11.94 


10.82 


7.9 


148.8 


1.75 




ESO239-IG002 


22 49 39.6 


-48 51 01 


0.043 


11.75 


10.88 


-0.1 


95.7 


1.13 




ESO244-G012 


01 18 08.6 


-44 27 40 


0.023 


11.32 


10.20 


5.3 


35.7 


0.42 




ESO264-G036 


10 43 07.0 


-46 12 43 


0.023 


11.24 


10.70 


1.5 


29.1 


0.34 




ESO286-IG019 


20 58 27.4 


-42 38 57 


0.043 


11.95 


11.13 


10.0 


151.6 


1.78 




ESO440-IG058 


12 06 53.0 


-31 57 08 


0.023 


11.33 


10.17 


99.0 


36.3 


0.43 




ESO507-G070 


13 02 51.3 


-23 55 10 


0.021 


11.44 


10.67 


6.8 


46.6 


0.55 




IC1623A/B 


01 07 46.3 


-17 30 32 


0.020 


11.63 


10.42 


6.0 


72.7 


0.86 


PSN 2011 


IC2545 


10 06 04.2 


-33 53 04 


0.034 


11.66 


10.58 


-0.1 


77.7 


0.92 




IC2810 


11 25 47.3 


+ 14 40 23 


0.034 


11.60 


10.61 


1.5 


68.0 


0.80 




IC4687/6 


18 13 38.6 


-57 43 36 


0.017 


11.45 


10.43 


2.9 


47.3 


0.56 


PSN 2010 


IRAS 12224-0624 


12 25 02.8 


-06 40 44 


0.026 


11.30 


9.83 


2.9 


34.0 


0.40 




IRAS 14378-3651 


14 40 57.8 


-37 04 25 


0.068 


12.13 


10.40 


5.1 


233.0 


2.74 




IRAS 16399-0937 


16 42 39.2 


-09 43 11 


0.027 


11.55 


10.39 


10.0 


60.3 


0.71 




IRAS 17207-0014 


17 23 21.4 


-00 17 00 


0.043 


12.42 


10.32 


-50. 


447.5 


5.27 




IRAS 18090+0130 


18 11 37.3 


+01 31 40 


0.029 


11.63 


10.54 


2.0 


71.7 


0.84 




MCG-02-01-051/2 


00 18 51.4 


-10 22 33 


0.027 


11.41 


10.59 


3.1 


44.4 


0.52 


2010 hp 


MCG-03-04-014 


01 10 08.5 


-16 51 14 


0.035 


11.59 


10.53 


-5.0 


66.1 


0.78 




NGC0034 


00 11 06.6 


-12 06 27 


0.020 


11.43 


10.34 


-1.0 


45.7 


0.54 




NGC0232 


00 42 46.5 


-23 33 31 


0.020 


11.51 


10.71 


1.1 


55.7 


0.66 


2006et 


NGC3110 


10 04 02.7 


-06 28 35 


0.017 


11.29 


10.94 


3.3 


33.0 


0.39 




NGC5010 


13 12 25.4 


-15 47 45 


0.021 


10.84 


9.79 


-1.0 


11.8 


0.14 




NGC5331 


13 52 16.6 


+02 06 08 


0.033 


11.60 


10.92 


3.0 


67.1 


0.79 




NGC6240 


16 52 58.6 


+02 24 03 


0.024 


11.81 


10.89 


-0.3 


108.4 


1.28 


2000bg, 2010gp 


NGC6926 


20 33 04.8 


-02 01 39 


0.020 


11.25 


11.38 


5.6 


30.4 


0.36 




NGC7130 


21 48 19.6 


-34 57 05 


0.016 


11.34 


10.75 


1.2 


37.1 


0.44 


2010bt 


NGC7592 


23 18 22.2 


-04 24 56 


0.024 


11.36 


10.51 


-1.0 


38.4 


0.45 




NGC7674 


23 27 56.9 


+08 46 46 


0.029 


11.37 


10.92 


1.1 


40.0 


0.47 


2000A, 2011ee, 2011hb 



the reference source list in general includes the bright galaxy nu- 
cleus which, being very bright, has a significant weight in the 
determination of the kernel. This may cause some problems be- 
cause if at one epoch a SN occurs very close to the galaxy nu- 
cleus it can be included in the convolution kernel and effectively 
cancelled in the difference image. We therefore modified the ISIS 
selection procedure to allow for exclusion of specific sources, in 
particular the galaxy nuclei, from the reference list. 

Despite the efforts in many cases the difference image shows 
significant spurious residuals in correspondence to the galaxy 
nuclear regions. The problem is most severe in case of im- 
ages with poor seeing (FWHM > 1") and/or reduced trans- 
parency. This is illustrated in Fig.[2]where we show two exam- 
ples of image difference one for a search image with poor seeing 
(FWHM = 1.5", left panel) and the other for a case with excel- 
lent seeing (FWHM = 0.4"). In both cases, the reference image 
was the same and had excellent seeing (FWHM = 0.4"). 
False detections due to residuals of the image subtraction were 
largely removed by the requirement that the candidate had to be 
visible at least in two consecutive epochs. 

2.4. Supernova discoveries and characterization 

During our monitoring campaign 6 transients were detected in 
at least two consecutive epochs separated by at least one month 
(finding charts are in Fig. |3}- Four of them were spectroscopi- 
cally confirmed as SNe (three SN-CC and one SN la) and we 
will argue in the following that also the other two transients, la- 
beled as probable SN (PSN), are likely SN CC (Tab. |3). SNe 
2010bt and 2010gp were discovered and announced before our 



W 




S3 









Fig. 2: Left panel: example of poor subtraction of images of 
NGC 7130 with large seeing differences (FWHM = 1.5", with 
FWHM - 0.4"' for the reference). Right panel: optimal subtrac- 
tion for two images with similar, good seeing (FWHM - 0.4"). 
In both panels the source in the lower left quadrant is SN 2010bt 
(cf. Fig.EJ. The FOV in both panel is about 2'x2'. 



detection by optical searches but have been independently re- 
discovered by us. 

The objects are listed in Tab. [3] along with the host galaxy 
name, distance modulus (computed from the galactocentric red- 
shift and the adopted cosmology), SN coordinates, offsets from 
the galaxy nucleus and projected linear distances from the 
galaxy nucleus. 

For all transients K-band magnitudes were measured through 
aperture photometry on the difference images and calibrated 
with respect to 2MASS stars in the field. Upper limits measured 
on pre-discovery images were also estimated. For all transients, 
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Table 2: The log of the observations with the epoch of observations (MJD), the seeing (FWHM in arcsec), and the minimum 
maximum magnitude limit for SN detection across the image (cf. Sec. 12.51) . 



JD 



CGCG011-076 

Seeing trfc 



ESO440-IG058 

JD Seeing m,'™ r 



IRAS17207-0014 



JD 




Seeing 


Urn 
max 


m'"" 

nun 


55035 





0.39 


19.0 


16.5 


55056 





0.37 


19.0 


16.5 


55090 


9 


0.56 


19.0 


15.3 


55342 


3 


0.78 


19.0 


14.2 


55398 


1 


1.16 


18.7 


16.7 


55415 


2 


0.49 


19.0 


14.8 


55446 


1 


1.05 


19.0 


16.0 


55675 


3 


0.42 


19.0 


18.0 


55739 


3 


1.05 


18.5 


17.0 


55775 





0.66 


19.0 


17.0 


55796 


1 


1.10 


18.7 


17.3 



JD 



NGC0232 

Seeing m,'™,. 



55189. 
55203. 
55219. 



0.41 
0.59 
0.65 



19.0 
19.0 
19.0 



13.5 
13.4 
13.1 



JD 



CGCG043-099 

Seeing m'Z x 



m 



lim 



54989.0 
55201.3 
55339.0 
55701.9 
55734.1 



0.36 
0.43 
0.56 
0.49 
1.8 



19.0 
19.0 
19.0 
19.0 
18.3 



15.5 
15.0 
17.3 
18.8 
15.5 



54989. 
55019. 
55238. 
55251. 
55284. 



JD 



0.49 
0.92 
0.98 
1.16 
0.42 
ESQ148- 
Seein: 



19.0 
19.0 
19.0 
19.0 
19.0 
1G002 



13.5 
13.0 
13.1 
12.5 
13.5 



ESO507-G070 



JD Seeing m 



lim 



m 



lim 



g m ;„ 



54989.1 
55016.1 
55226.2 
55251.3 
55284.0 



0.36 
0.53 
0.36 
0.53 
0.47 



19.0 
19.0 
19.0 
19.0 
19.0 



16.5 
16.7 
16.5 
16.5 
16.7 



IRAS18090+0130 



54990. 
55028. 
55048. 
55098. 
55365. 
55375. 
55393. 
55702. 
55719. 
55743. 



0.68 
0.69 
0.42 
0.57 
0.56 
1.02 
0.99 
0.71 
1.47 
0.59 



19.0 
19.0 
19.0 
18.7 
18.8 
19.0 
19.0 
19.0 
18.2 
18.8 



15.0 
15.0 
15.0 
15.0 
15.3 
14.0 
15.0 
16.8 
15.0 
16.8 



JD 



55056.3 
55071.3 
55390.4 
55415.3 
55444.3 
55762.3 
55796.3 
55812.1 



IC1623A 

Seeing m,'™ r 



m 



0.41 
0.67 
0.54 
0.37 
0.51 
0.89 
0.79 
0.63 



19.0 
19.0 
19.0 
19.0 
19.0 
19.0 
19.0 
18.8 



15.5 
15.5 
17.3 
16.6 
17.3 
13.5 
15.3 
14.8 



JD 



ES0239 

Seein: 



IG002 



m : 



JD 



IC2545 

Seeing m 



lim 
mux 



111 



JD 




Seeing 


lim 
max 


m'"" 

nun 


55019 


1 


1.19 


19.0 


16.0 


55038 


1 


0.51 


19.0 


16.7 


55056 





0.65 


19.0 


16.2 


55340 


3 


1.00 


18.7 


13.2 


55340 


3 


1.00 


18.7 


13.2 


55366 


3 


0.58 


19.0 


15.5 


55398 


1 


1.20 


18.2 


14.8 


55415 


2 


0.41 


18.8 


16.8 


55672 


3 


0.39 


19.0 


15.7 


55742 


2 


0.47 


19.0 


15.3 


55796 


1 


0.89 


18.5 


15.3 



55019. 
55028. 
55056. 
55342. 
55365. 
55398. 
55415. 
55675. 
55702. 
55719. 



0.72 
0.57 
0.53 
0.63 
0.69 
1.16 
0.37 
0.50 
0.71 
1.58 



19.0 
19.0 
19.0 
19.0 
19.0 
18.7 
19.0 
19.0 
18.8 
18.2 



14.7 
14.5 
14.5 
14.5 
15.5 
14.5 
14.7 
15.0 
14.5 
14.8 



55137.3 0.95 

55172.3 0.48 

55662.0 0.56 

55665.1 0.57 



19.0 
19.0 
19.0 
19.0 



15.5 
15.5 
16.0 
15.5 



MCG-02-01-051 



JD 



IC28HI 

Seeing m'Z x 



m 



lilll 



55202.3 0.50 

55226.2 0.51 

55238.3 0.46 



19.0 
19.0 
19.0 



14.0 
13.8 
14.8 



JD 



ES0244 

Seein: 



G012 

lim 



JD 



g m; 



m 



lim 



IC4687 

Seeing m,'™ r 



JD 


Seeing 


m'"" 

"max 


m'"" 

nun 


55037.4 


0.56 


19.0 


18.5 


55049.4 


0.47 


19.0 


18.8 


55068.4 


0.47 


19.0 


18.2 


55398.2 


0.92 


19.0 


18.8 


55415.3 


0.45 


19.0 


19.0 


55449.3 


0.47 


18.7 


18.5 


55704.3 


1.20 


18.7 


18.3 


55739.2 


1.05 


19.0 


18.7 



MCG-03-04-014 



55028 
55056 
55071 
55375 
55393 
55402 
55705 
55757 



1.11 
0.39 
0.63 
1.16 
1.20 
0.82 
1.26 
0.84 



19.0 
19.0 
19.0 
18.0 
18.2 
18.8 
18.5 
19.0 



13.0 
14.0 
13.5 
13.7 
14.0 
13.5 
15.5 
16.2 



54989.1 
55056.0 
55337.1 
55366.3 
55375.0 
55393.2 
55663.3 



0.37 
0.65 
0.55 
0.58 
1.21 
1.21 
0.66 



19.0 
19.0 
19.0 
19.0 
19.0 
19.0 
19.0 



16.0 
15.8 
16.2 
16.2 
16.2 
15.3 
16.2 



JD Seeing m;™^ mj 



1RAS12224 

JD Seeing 



TJ52T" 

lim 



m 



55056 
55106 
55443 
55762 
55796 
55812 
55402 
55415 



0.48 
0.74 
0.84 
1.05 
0.63 
0.53 
1.16 
0.39 



19.0 
19.0 
19.0 
19.0 
19.0 
19.0 
19.0 
19.0 



JD 



ES0264 

Seein: 



G036 

Jim 



g m m 



55142. 
55172. 
55184. 
55339. 
55665. 
55701. 
55734. 



0.98 
0.63 
0.45 
0.59 
0.55 
0.50 
1.80 



18.0 
18.3 
18.7 
18.3 
19.0 
19.0 
18.8 



14.0 
14.5 
14.0 
13.5 
16.8 
17.5 
16.7 



54989.0 0.51 
55019.0 0.84 
55226.2 0.47 



19.0 
19.0 
19.0 



16.5 
16.0 
16.0 



NGC0034 



IRAS 14378 

JD Seeing 



3651 

m lim 



54990.1 0.53 
55072.9 0.36 

1RAS16399 

JD Seeing 



19.0 
19.0 
TJ93T" 



13.5 
14.5 



Jim 



JD 



ES0286 

Seein: 



IG019 

liiii 



g m ;„ 



Jim 



54988 
55019 
55038 
55056 
55341 
55365 
55382 
55398 
55700 
55719 
55735 



0.34 
0.62 
0.45 
0.45 
0.74 
0.49 
0.63 
0.84 
0.65 
1.34 
1.60 



19.0 
19.0 
19.0 
19.0 
19.0 
19.0 
19.0 
18.7 
19.0 
18.5 
19.0 



15.0 
15.5 
15.0 
15.5 
14.5 
15.2 
15.0 
15.0 
16.8 
15.2 
14.5 



55014.1 
55035.0 
55057.0 
55091.0 
55342.3 
55398.0 
55416.1 
55446.0 
55665.2 
55710.0 
55771.1 
55785.0 
55795.0 



0.49 
0.42 
0.38 
0.59 
0.80 
1.20 
0.40 
1.00 
0.60 
0.63 
0.37 
0.89 
0.87 



19.0 
19.0 
19.0 
19.0 
19.0 
18.7 
19.0 
19.0 
19.0 
18.7 
18.3 
18.8 
18.9 



16.0 
16.7 
16.3 
16.8 
15.7 
15.7 
16.5 
16.3 
17.0 
16.3 
15.7 
16.2 
16.2 



NGC7674 



JD 



NGC0232 

Seeing m', 



lim 
mux 



55090. 
55382. 
55402. 



0.37 
0.69 
0.95 



19.0 
19.0 
19.0 



-.lim 



17.0 
16.2 
14.8 
15.8 
14.8 
15.0 
14.5 
14.8 



JD 


Seeing 


Urn 

IIIO.X 


m'"" 

nun 


55049.3 


0.47 


18.5 


14.3 


55068.3 


0.53 


18.8 


14.0 


55398.2 


1.37 


18.2 


12.7 


55449.3 


0.40 


19.0 


13.5 


55704.4 


0.79 


19.0 


14.0 


55739.2 


1.00 


19.0 


15.3 


55771.3 


0.59 


18.8 


14.2 



JD 


Seeing 


lim 

max 


m 

nun 


55068.1 


0.47 


19.0 


14.5 


55398.1 


1.24 


18.0 


12.7 


55415.3 


0.42 


19.0 


13.8 


55446.1 


0.84 


19.0 


14.0 


55739.3 


0.68 


19.0 


15.7 


55796.3 


0.74 


19.0 


15.3 



111 



15.2 
14.7 
16.0 



55415.3 
55704.3 
55762.3 



0.34 
1.02 
1.05 



19.0 
19.0 
19.0 



16.0 
14.0 
14.2 



JD 



NGC3110 

Seeing m',Z x 



m: 



lim 



55184.3 
55202.3 
55291.0 
55338.9 
55662.0 
55708.0 
55189.3 



0.47 
0.35 
0.45 
0.56 
0.68 
0.42 
0.33 



19.0 
19.0 
19.0 
19.0 
19.0 
19.0 
19.0 



16.0 
16.8 
15.8 
16.7 
15.0 
15.5 
16.2 



JD 



55014.9 
55035.0 
55243.3 
55284.0 
55402.0 
55665.1 
55710.0 



NGC5010 

Seeing mj,™. 



m 



0.50 
0.58 
1.02 
0.44 
0.99 
0.48 
0.65 



19.0 
19.0 
18.3 
19.0 
19.0 
19.0 
19.0 



16.0 
14.5 
14.2 
17.7 
14.5 
16.3 
16.5 



JD 



NGC5331 

Seeing 



Jim 



mlim,,. 



55016.1 
55251.3 
55263.3 
55665.2 
55708.0 



0.49 
0.95 
0.48 
0.49 
0.42 



19.0 
18.7 
19.0 
18.8 
17.7 



14.0 
14.8 
15.0 
16.3 
16.5 



JD 



NGC6240 

Seeing m^' x 



m 



lim 



55019.1 
55038.1 
55057.0 
55342.3 
55398.0 
55416.1 
55455.9 
55771.1 
55785.1 
55796.0 



0.54 
0.53 
0.45 
0.82 
1.00 
0.36 
0.55 
0.47 
1.03 
1.26 



19.0 
19.0 
19.0 
18.3 
17.7 
18.2 
18.7 
19.0 
19.0 
19.0 



14.7 
17.5 
15.8 
12.3 
12.3 
11.0 
13.0 
15.0 
14.7 
13.5 



JD 



NGC6926 

Seeing 



m' 



lim 



55068.1 
55366.3 
55398.1 
55415.2 
55672.3 
55701.3 
55743.1 



0.44 
0.43 
1.26 
0.34 
0.39 
0.74 
0.45 



19.0 
19.0 
18.5 
19.0 
19.0 
18.7 
19.0 



14.8 
15.0 
14.0 
16.3 
16.7 
15.5 
17.7 



JD 



NGC7130 

Seeing 



55038.1 
55341.3 
55366.3 
55382.2 
55675.3 
55719.2 
55735.2 



0.47 
0.76 
0.44 
0.57 
0.40 
1.37 
1.0 



19.0 
19.0 
19.0 
19.0 
18.5 
18.2 
19.0 



16.0 
14.7 
16.3 
16.7 
13.7 
15.3 
16.0 



JD 



55398.2 
55415.2 
55444.2 
55704.3 
55771.3 
55796.3 



NGC7592 

Seeing 



m 



1.0 
0.37 
0.40 
0.95 
0.47 
0.65 



17.5 
18.7 
19.0 
19.0 
19.0 
19.0 



13.8 
14.2 
15.0 
13.5 
15.5 
13.5 
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Table 3: Information for the detected SNe. 





SN2010bt 


SN2010gp 


SN2010hp 


SN201 lee 


PSN2010 


PSN201 1 


Host galaxy 


NGC7130 


NGC6240 


MCG-02-01-52 


NGC7674 


IC 4687 


IC 1623 A 


Distance modulus 


34.08 


35.07 


35.24 


35.39 


34.21 


34.54 


R.A. (SN) 


21:48:20.22 


16:52:57.39 


0:18:50.01 


23:27:57.34 


18:13:40.213 


01:07:46.229 


Dec. (SN) 


34:57:16.5 


2:23:16.4 


-10:21:40.6 


+08:46:38.10 


-57:43:28.00 


-17:30:29.48 


Offset ["] 


9E 14S 


22W 47 S 


2.6W 2.7N 


9.3E 6.3S 


2.6E 2.8N 


4E 7S 


r [Kpc] 


5.3 


25.9 


2.1 


8.9 


6.5 


3.2 


SN type 


Iln 


la 


IIP 


Ic 


IIP* 


Ic* 


R max [mag] 


i j.y + U. 1 


i A n _i_ n i 
10. U ± U.l 


i h a _i_ n 1 
1 l.b ± U.Z 


1 /. j ± U.j 




1 A o _i_ n c 


K max [mag] 


16.1+0.2 


16.9 ±0.1 


17.2 ±0.1 


18.0 ±0.5 


17.9 ±0.2 


17.0 ±0.3 


AB G [mag] 


0.44 


0.33 


0.16 


0.25 


0.43 


0.07 


AB H [mag] 


1.7 ±0.5 


0.2 ±0.1 


0.5 ±0.3 


-0.0 


0-8 


0.5 ±0.5 


abs M R 


-19.2 


-19.2 


-17.9 


-18.1 




-17.9 


abs M K 


-18.1 


-18.2 


-18.0 


-17.7 


-16.3/ -17.3 


-17.5 


MJD max R 


55303 ± 5 


55405 ± 3 


55399 ± 5 


55760 ± 5 


55270 ± 20 


55725 ± 10 


Refs 


CBET 2250, 2252 


CBET 2388, 2390 


CBET 2246, 2249 


CBET 2773 







# Photometric classification. 



but PSN2010 in IC4687, we obtained some follow-up imaging 
in the optical or near-infrared domains. These observations were 
reduced using standard procedures in IRAF. When a reference 
image was not available, the SN magnitude was measured using 
the PSF fitting technique. Optical band magnitudes were cali- 
brated with respect to Landolt's standard fields. Our photometry 
for the six transients is reported in Tab - ] For the two transients 
with no spectroscopic confirmation, the photometry will be used 
to assess their nature. 

Spectroscopic observations were obtained for four candi- 
dates: epoch, spectral range and instruments are reported in 
Tab - ] Data were reduced using standard procedure in IRAF 
but for the X-Shooter spectra which wer e reduced using ver- 
sion 1.0.0 of the ESO X-shooter pipeline dGoldoni et al.l l2006) 
with the calibration frames (biases, darks, arc lamps, and flat 
fields) taken during daytime. The extracted spectra, after wave- 
length and flux calibration, were compared with a library of tem- 
plate spectra using the GELATO SN spectra co mparison tool 
(https://gelato.tng.iac.eS/. lHarutviinvan~e t al. 2008). The best fit 
template SN, the SN type and phase are reported in Tab - ] 
Spectroscopic classification for PSN2010 in IC4687 was at- 
tempted, but the observed spectrum resulted too noisy for a safe 
classification. The table includes t he result of the spe ctroscopic 
observations of SN 2010gp from lFolatelli et all (120101) . 

We have used the available photometry and spectroscopy to 
put some constraints to the amount of extinction suffered by the 
SNe. Hereafter we will describe in some details the sparse infor- 
mation available for each transient. 



SN 2010bt was discovered on 2010 April 17.10 UT by Monard 
(I2OIOI) . A spectrum taken on April 18.39 UT dTuratto et al. 
2010) shows strong r esemblance to seve ral type-IIn SNe, in 
particular SN 1996L dBenettiet all 19991) shortly after explo- 
sion (Fig. - }. A broad H a component is present indicating an 
expansion velocity of about 3500 kms~' (half width at zero 
intensity). SN 2010bt was independently re-discovered by 
us on May 25 and was observed in other 2 epochs. The ob- 
ject was not visible on a HAWK-I image taken in 2009 July 
26 (limit K - 19.0 mag). While the analysis of the light 
curve and spectral evolution of th is SN will be presented 
elsewhere (El ias-Rosa et al.l 1201 3b . a preliminary analysis 
shows that to match the color of SN 2010bt to that of the 
type Iln SN 1998S requires a significant amount of extinc- 
tion (Ab = 1.7 ± 0.5 mag). This extinction improves also the 
matching of the spectra of SN 2010bt to SN 1996L. 



SN 2010gp wa s discovered on 2010 July 14.10 UT by 
iMaza et alJ <l2010h with the 0.4 1-m PROMPT ! telescope lo- 
cated at Cerro Tololo. iFolatelli et al.l (1201 Ol) reported the 
spectroscopic classification as a type-la SN around maxi- 
mum light and with high expansion velocity of the ejecta. 
SN 2010gp was re-discovered independently by us on July 
21 and was observed in other 2 epochs. The object was not 
visible on a HAWK-I image taken on 2010 May 26 (limit 
K = 18.5 mag). The available colors compared to that of 
standard SN la suggest a small reddening inside the host 
galaxy, A B = 0.2 mag. 

SN 2010hp was dis covered on a HAWK-I ima ge taken on 2010 
July 21.3 UT dMiluzio & Cappellardl20T~l) . The object was 
not detected on 2009 Aug. 25 (limit K = 19.0 mag). Based 
on a spectrum taken on 2010 Sept 8, the SN was classified 
as a type II event m ore than 30 days past maximum light 
(Mari on et al.l l2010). We obtained a follow-up spectrum on 
2010 Sept 15 with EFOSC2 at the NTT showing the typical 
features of type II SN, that is H Balmer lines, Nal D dou- 
blet, the NIR Call emission triplet along with a number of 
Fell lines in the blue region. The GELATO spectral compar- 
ison tool found a best match with the type IIP SN 1999em 
(Elmhamdi et a~ l2003l) at about +60 days, adopting a red- 
dening of about = 0.5 mag (Fig - }. This is fully consis- 
tent with the color curves comparison between the two SNe. 

SN 2011ee was discov ered on 2011 June 27.3 UT 
(Mil uzio et al.l l2fjj"~) . The object was not detected on a 
K-band image taken on 2010 September 7 (K > 19.0 mag). 
An optical spectrum was obtained with X-Shooter at the 
VLT on 2011 July 17.3 UT showing that the transient 
is a type Ic SN. The GELATO co de finds a best match 
with SN 2007gr dHunteret al.l 120091) at maximum (Fig. @). 
Because the classification of type Ic SNe near maximum 
is sometimes ambiguous we obtained a followup spectrum 
about two months later (2011 Sept 20) using OSIRIS at 
the GTC. Again the spectrum is very similar to SN 2007gr 
at corresponding phase (Fig. 0] bottom panel). The color 
comparison and the spectral match are consistent with a 
negligible host galaxy extinction. 

PSN2010 in IC 4687 was discovered on 2010 May 21.3 UT in 
the northern component of a galaxy triplet that include also 
IC 4686 and, 1 arcmin to the south of IC 4687, IC 4689. 
IC 4687 has a chaotic structure formed byf stars, gas and 
dust and a large curly tail. The transient was not detected on 
a K-band image taken on 2009 Aug. 8.1 ( K > 19.0 mag). 
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Fig. 3: K-band finding charts for the SNe of our list. The inserts show the transients as they appear in the difference image. 
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Table 4: Transient photometry. Estimated errors are given in parentheses. 



MJD 


B 


V 


R I 


J 


H 


K 


Instr. 








SN 2010bt 










55341.4 












16.9 (.1) 


H 


JJ J 00. 4 












i o no { 1 \ 
lo.Us (.1) 


rl 


55382.3 












18.68 (.1) 


H 








SN 2010gp 










55342.4 












> 18.5 


H 


55398.1 












17.4 (.1) 


H 


55416.1 












16.85 (.1) 


H 


55457.0 












17.9 (.1) 


H 


55455.0 


19.72 (.05) 


18.43 (.03) 


18.11 (.03) 18.15 (.03) 








E 


55398.3 






2010hp 






17.25 (.05) 


H 


55415.4 












17.17 (.05) 


H 


55449.4 












17.20 (.10) 


H 


55454.3 


19.61 (.03) 


18.70 (.03) 


18.28 (.03) 17.89 (.03) 








E 


55456.1 








17.24 (.1) 


16.93 (.1) 


17.17 (.1) 


S 


55498.1 








18.48 (.1) 


18.51 (.1) 


17.99 (.1) 
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H = HAWK-I@VLT, E = EFOSC2@NTT, S = SOFI@NTT, L = RATCam@Liverpool, D = Dolores® TNG 



Table 5: Log of spectroscopic observations. 



SN 


MJD 


range(nm) 


res. (nm) 


Instrument 


best fit 


type Phase 


2010bt 


55304.4 


350-1000 


1.4 


EFOSC/NTT 


2005gj 


Iln max 


2010hp 


55454.2 


350-1000 


1.4 


EFOSC/NTT 


1999em 


la +60d 


2010 IC4687 


55352.2 


340-1000 


0.2 


XSHOOTER/VLT 




S/N too low 


201 lee 


55824.6 




2.5 


OSIRIS/GTC 


2007gr 


Ic +60 


201 lee 


55759.3 


340-1000 


0.2 


XSHOOTER/VLT 


19941 


Ic max 


2010gp 


55401.0 


354-886 


Folatelli et al. (2010) 


2002bo 


la just before max 



We obtained an optical/infrared spectrum with X-Shooter at 
VLT on 2010 June 5. However, because of its very low S/N, 
we could not derive a convincing classification and there- 
fore we had to rely on the K-band photometry. Comparing 
the K band absolute light curve of PSN2010 (A B (host) = 0) 
with template light curves of different SN types we found a 
good match with SN 2005cs a pro totype of under-luminous 
type IIP SN (iPastorello et alJl20 09). assuming that the detec- 
tion of PSN2010 was 2 months after the explosion. However, 
lacking color measurements, we could not constraint the ex- 



tinction and indeed, assuming a high extinction Ag ~ 8 mag, 
we found an alternative good match with the light curve 
of SN 1999em (Fig. Q. Intermediate values may also be 
adopted by fitting other SN II. 
PSN2011 in IC 1623 was discovered on 201 1 July 21.4 UT in 
the western component of a galaxy pair. The object was not 
detected on a K-band image taken on 2010 Sept. 5 (K > 19.0 
mag). Unfortunately, due to bad weather in the scheduled 
nights, we could not obtain a spectroscopic observation of 
the transients. We have to rely on three epochs of photome- 
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(a) The spectrum of 2010bt, dereddened by A B = 1.7 mag 
(see text) is compared with that of the type Iln SN 1996L. 
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(b) The spectrum of 2010hp, dereddened by A B = 0,5 mag 
(see text) is compared with that of the type IIP SN 1999em. 
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(c) The spectrum of 201 lee is compared to that of the 
SN Ic 2007gr at the maximum (top panel) and 60 days 
after the maximum. 



Fig. 5: K band absolute light curve of the PSN2010 (black dots), 
compared with those of template SNe. The empty circles show 
the same data but assuming an extinction Ag = 8 mag. 

try, in K complemented by two epochs in the optical R and I 
bands. A simultaneous comparison of the absolute observed 
luminosity with template SNe give a best fit with the SN Ic 
2007 gr one month after maximum (Fig. |6}. Assuming this 
classification, from the colors we can constrain the extinc- 
tion to be Ab = 0.5 + 0.5 mag. However, we have to admit 
that, within the errors, the photometry of PSN2011 can be 
consistent also with a type IIP at about 3 months after explo- 
sion. 

To recap, during the search we discovered 6 SNe. Four re- 
ceived a spectroscopic classifications: one as a type la and three 
as core collapse events, a type Iln, a type IIP and a type Ic. For 
the other two, based on the sparse photometry, we argue that 
most likely they are core collapse SNe, with a best fit as type IIP 
and type Ic, respectively. 

2.5. Search detection limit 

In order to derive the SN rate from the number of detected events 
it is crucial to obtain an accurate estimate of the magnitude de- 
tection limit for each of the search images and for different loca- 
tions in the images. As it has been shown in Fig. [2 the detection 
efficiency is influenced by the sky conditions at the time of ob- 
servations (namely seeing and transparency) and by the transient 
position inside the host galaxy. 

The magnitude limit for SN detection has been estimated 
through artificial star experiments. The procedure we adopted 
was the following: 

1 . fake SNe of different magnitudes are simulated with the PSF 
derived from isolated field stars; 

2. the image is segmented in a number of intensity contour lev- 
els. We took denser contours in the nuclear regions because 
the magnitude limit changes rapidly with background inten- 
sity; 

3. one fake SN of specific magnitude is randomly placed inside 
a chosen intensity contour; 



Fig. 4: Spectra of SN 2010bt, 2010hp, and 201 lee are shown 
along with the best fitting templates. 
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Fig. 6: K and R-I colours of PSN 201 1, compared with those of 
template SNe. For this plot we adopted Ab — 0.5 mag, which 
improves the fit of the R band photometry (AR = 0.3 mag). 

4. the image with the fake SN is processed through the image 
difference and transient detection pipeline; 

5. if the fake SN results in a detected transient, the experiment 
is repeated with a fainter artificial star until we have a null 
detection. The fainter magnitude for which the fake SN is de- 
tected defines the magnitude limit for the given background 
intensity level; 

6. Steps 2 to 5 were repeated for each contour level three times 
to enhance the statistical significance of the results. The av- 
erage value for each contour level has been adopted as the 
magnitude discovery limit for the given background inten- 
sity. 

To illustrate the results, a plot of the magnitude limit ver- 
sus background counts for four observations of the galaxy 
NGC 7 130 is shown in Fig. Q Each epoch is labelled with the im- 
age seeing, while the errorbar shows the range of limiting magni- 
tudes for the three experiments. The top x-axis shows the linear 
distance in Kpc from the galaxy center. 

It can be seen that, as expected, the magnitude limit is lower 
in the nuclear regions which, for a typical galaxy, correspond to 
1.5-2.0 kpc. Epochs with different seeing have similar magnitude 
limits in the galaxy outskirts (typically K ~ 19 mag), while in 
the nuclear region when seeing is poorer the magnitude limit 
is brighter (in the worst case even 5-6 mag brighter than in the 
galaxy outskirts). 

3. SN search simulation 

To evaluate the significance of the detected events we elabo- 
rated a simulation tool that returns the number and properties 
of expected events based on specific features of our SN search, a 
number of parameters describing our current knowledge of SBs 
and SN properties. The tool uses a MonteCarlo approach which 



Kpc 

>5 2 1 0.3 
19-5 i 1 1 1 1 1 1 1 r 



E 




2000 4000 6000 8000 10000 12000 14000 



counts 

Fig. 7: Magnitude limit vs. background intensity for four differ- 
ent observations of NGC 7130. The points are labeled with the 
respective image seeing. Errorbars show the range of the magni- 
tude limits from the three different experiments. 

simulates the stochastic nature of SN explosions. By collecting 
a number of MonteCarlo experiments with the same input pa- 
rameters, we can test whether the observed events are within the 
expected distribution. On the other hand by varying some of the 
input parameters, we can test the influence of specific assump- 
tions. 

3.1. The simulation tool 

Our MonteCarlo (MC) simulation tool is built in a Python en- 
vironment and makes use, for the different inputs, of standard 
values taken from the literature. For those that are more contro- 
versial, we will give references with some discussion. The basic 
ingredients of the simulation are: 

- relevant data for the selected SBs, namely: redshift, galac- 
tic extinction, infrared fluxes from IRAS catalogues at 25, 
60 and 100//, B magnitude corrected for internal extinction^, 
Hubble morphological type. These data were retrieved from 
NED; 

- information describing the SN properties for each of the SN 
types considered here: SNe la and core collapse events, in- 
cluding SNe IIP, ILL, Iln and Ib/c. K band template light 
curves were constructed star ting from B template light 
curves (Cappe llaro et al.ll 19971) and B — K k-corr ections for 
the given galaxy redshift dBotticella et al.l 12008). We note 
that the results of this procedure are in close agreemen t 
with the K template light curves of Mattila & Meikld <2001l) . 
For the SN luminosity f unctions we adopted as reference 
those of iLi et all (1201 lbl) . but we also tested for the possi- 
ble presence of a significant population of faint core collapse 
which may be s uggested by the analysis of very nearby SNe 
dHoriuchi et al.ll201lh . From the LOSS project we adopted 
also th e relative rates for the different SN types (ILi et al.l 
l20lT3> : 

- details of the search campaign: log of observations, magni- 
tude detection limit for each observation as a function of the 
host galaxy background intensity; 

- the number of SNe expected from a given star formation 
episode. For core collapse SNe, this is determined only by 

4 The galaxy in ternal extinction was retrieved from HyperLeda 
iPatur el et al. 20Q3) 
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the adopted mass range of the progenitors and IMF slope. 
In fact, for our purposes, we can neglect the very short time 
delay from CC progenitor formation to explosion. For type 
la SNe we need to consider the realization factor, that is the 
fraction of events in the proper mass range which occurs in 
suitable close binary systems and the delay time distribution 
(Sect [3X3.; 

- the depth and distribution of the extinction by dust inside the 
parent galaxies (cf. Sect. 13.1. 3k 

- the star formation spatial distribution in the parent galaxies 
(cf. Sect. EH. 

Hereafter, we discuss our assumptions about the parameters 
of the simulation. 

3.1.1. From IRAS measurements to Star Formation Rate 

The SFR in SBs can be estimated on the basis of the galaxy 
total infrared luminosity (Ltir) under the assumption that dust 
re-radiates a major fraction of the UV luminosity, and after cali- 
bration with stellar synthesis models. In turn the TIR luminosity 
can be estimated from FIR flux measurements. 

iHelou et"aT] ([1988) provided a prescription for deriving the 
FIR emission from IRAS measurements: 

FIR = 1.26 x 1(T 14 [2.58/ v (60//m) + / v (100yuw)] 

where FIR is in W m~ 2 and f v are in Jansky. FIR fluxes 
are co nverted into TIR fluxes by using the relation o f lDale et al. 
J2001 

1 2 3 4 

log = fly + a\x + CI2X + cit,x + CI4X 

FIR 

where x = log /ff^ and [a(z = 0)] [0.2378, -0.0282, 
0.7281,0.6208,0.9118] 

TIR fluxes are converted in luminosities using the adopted 
distances: 

L TIR = 4nD 2 TIR 

Final ly, the relat i on bet ween the SFR (if/), and Ltir was de- 
rived by iKennicutl (Il998l) from SB galaxy spectral synthesis 
model adopting 10-100 Myr continuous bursts and a Salpeter 
IMF as: 

ft _ Ltir Ltir 

[Moyr- 1 ] ~ 2.2 x 10 43 [ergs- 1 ] ~ 5.8 x 10 9 [L ] 

3.1.2. SNRand SFR 

In general, the rate of SNe expected at a specific time, hsN(t), 
for a stellar population depends on the star formation history, 
the number of SNe per unit mass from one stellar genera- 
tion (labelled as SN productivity) and the distribution of delay 
time from star formation to explosion for the s pecific SN type. 
Following the notation of lGreggiol (I20051 l2010h : 

hsN(t)= f ^{t-r)k SN f SN (T)dr (1) 
Jo 

where ifr(t) is the star formation rate, /sn is the distribution 
of the delay times r and ksN is the supernova productivity. The 



equation shows that at a fixed epoch t since the beginning of star 
formation, the rate of SNe is obtained by adding the contribution 
of all past stellar generations, each of them weighted with the 
SFR at the appropriate time. 

Core Collapse SNe 

For core collapse SNe the delay time from star formation to ex- 
plosion (2.5 Myr for 120 M G stars up to 40 Myr for 8 M stars) 
is short compared wit h the typical SBs duration (200 - 400 Myr, 
McO uinn et alj|2009l) . Assuming that the SFR in the SB was 
constant during the past 40 Myr, the expected CC SN rate, hcc, 
is proportional to the current SFR: 

«cc = k C c X ^ 

The supernova productivity kcc is derived by integrating the 
IMF, <p(m), and assuming a CC progenitor mass (Mcc) range: 

Jug ^ m)dm 

kcc ~ r M » M w 

J M m<p(m)am 

where M cc and M^ c are respectively the lower and upper 
mass limits for SN CC progenitors and Ml, Mjj are the lower and 
upper stellar mass limit. To be consistent with the Kennicutt's 
SFR calibration we adopted a Salpeter IMF, that is: 



<f>(m) oc m~ a with a = 2.35 and 0.1 M Q < M < 100 M Q 

Assuming 8 < Mcc < 50 M0 for the CC progenitor mass 
range, kcc - 0.007 M _I . This number changes significantly if 
we adopt a different IMF, e.g. kcc = 0.01 1 f or a Kroupa IMF o r 
k C c = 0.039 for an extreme Starburst IMF dDwek et alJl201 lh . 
We soon note however that, because the IMF enters also in the 
conversion from Ltir to iff, the expected rate of SN events is 
almost independent on the selected IMF (cf. Sect. 14.11 ) provided 
the choice is consistent. 

More important is the assumption on the mass range for 
CC progenitors which is not well constrained. Actually, while 
changing the upper limit of the progenitor mass from 40 to 
100 M0 makes a modest 10% increase in the CC SN produc- 
tivity, the lower mass limit is crucial, with kcc decreasing by 
30% if w e adopted M^ c = IOMq instead of the favored value 
of8M^ (ISmarttll2009l) . 

SN la 

Estimating the expected rate of SN la is complicated be- 
cause the delay time distribution fi„, while still uncertain, cer- 
tainly ranges from short to very long time. In particular it has 
been suggested that SN la can be divided into two classes, 
one with a short delay time whose rate scales with the cur- 
rent SFR (also called prompt), and a second with a long de- 
lay time (tardy), whose rate scales with the average of the 
SFR along the entire galactic evolution (Scannapieco & Bildsten 
200 51 iMannucci et alJl2 006). Whil e stellar evo lution arguments 
(Greggio 2010, 2005; Greggio & Renzini 1983) and more recent 
data (Mao z et al.ll2012tlTotani et al.ll2008l) suggest a continuous 
distribution of the delay time instead of two distinct classes, the 
schematization is still a fair approximation that help in simplify- 
ing the problem of predicting the expected SN la rate in SBs. 
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In general, for a galaxy of the local Universe, ~ 13 Gyr after 
the beginning of SFR, we can identify th e contribution of the 
two components as follows dGregg io 2010): 

A/ B (13) = k la x J<Ac JT f la (r) dT + i/f P f f la {T) ofrj (2) 

where ij/c and tffp are the average SFR over, respectively, the 
last 0.1 Gyr (current SFR) and from 0.1 to 13 Gyr ago (past 
SFR). The SN productivity kj a is the product of the number of 
stars per unit mass in the adopted progenitor mass range (0.021 
for a Salpeter IMF and a mass range 3M0 < M < 8 Mq) and the 
realization fraction, the actual fraction of systems which make 
a succes sful explosion (~ 5% acc ording to the most recent es- 
timate) (Maoz & Mannuccill2012h . We assume that SF history 
in SBs can be described schematically with two components: 
a constant SFR during the galaxy evolution which created the 
galaxy stellar mass, and an on-going episode of intense SFR 
which is the source of the strong TIR emission. Neglecting the 
contribution of the ongoing SB to the galaxy stellar mass, we 
can approximate t//p M/13 x 10 9 , and write Eq.|2]as follows: 

«/a(13) - kja [tycFj a + J^yF'la) 

where F p Ia =< f la > xO.l and F\ a =< f Ia > Xl3 - 1 - F* a 
are the relative fraction of prompt and tardy events derived by 
integrating the delay time distribution in the relevant time range. 
In our approximation tffc can be derived from the observed Ljip 
and the galaxy mass from the K magnitude and B - K colors (cf. 
iMannucci etal]|2005l) . 

The relative contribution of the two SN la components has 
been a very debated issue in the last few years, ranging from 
F p Ia ~ 50% (IMannucci et alJl2006h to F p . ~ 10% from standard 
stellar evolution scenarios (lGreggioll2010l) . In our simulation we 
adopted as reference an intermediate value, F? ~ 30%. 

3.1.3. Extinction 

Dust extinction in S Bs is very high, esp ecially in the nuclear 
regions. For instance. IShiova et al.l (l200ll) found that fitting the 
spectral energy distribution of the nuclear region of Arp 220 re- 
qu ires a visua l extin ction Ay > 30 mag. Actually, according 
to lEngel et ail (1201 lh . "over most of the disk the near-infrared 
obscuration is moderate, but increases dramatically in the cen- 
tral tens of parsecs of each nucleus". Similar high extinction. 
Ay ~ 20, was found for the SB region of Zw 096 (Inami eTaTI 
l2010h . 

As a first order approximation, for our simulation we as- 
sumed that the extinction has the same distribution of the SF (see 
next section) with a maximum value Ay = 30 mag correspond- 
ing to the SFR peak and scaled linearly in the other regions. 
While this is a crude approximation, it turns out that the actual 
choice of extinction correction has little impact for our simu- 
lation. In the nuclear, high extinction regions the SN detection 
is limited by the reduced performance of the image subtraction 
algorithm in these high surface brightness regions. At the same 
time, our IR search is largely insensitive to variation in the (mod- 
erate) extinction of the outer galaxy regions. 

For the wavelength dependence o f extinction we ado pted the 
Calzetti's law with R v = 4.05 + 0.8 dCalzetti et alJl2000l) . 



3.1.4. Star Formation Distribution 

The spatial distribution of the SFR is a key ingredient of the 
simulation. This is because we expect that SNe occur more fre- 
quently in the high SF regions where, on the other hand, our de- 
tection efficiency is lower. In principle, the FIR emission which 
is used to estimate the SFR would also be a good tracer of its spa- 
tial distribution. However, it turned out that the available MIR 
imaging for the galaxies of our sample (mainly obtained with 
the Spitzer observatory) do not have enough spatial resolution 
for mapping the compact SB structures. 

Selected K-band images from our survey can have excellent 
resolution but, as is well-known, the near IR emission better 
traces the old star population, that is the galaxy mass distribu- 
tion more than the SFR distribution. Therefore for an estimate 
of the SFR concentration, we are forced to an indirect, statistical 
approach. 

Ou r starting point is the SB classification by lHattori et al.l 
(2004), who derived a correlation between the global SBs prop- 
erties, such as FIR colors, and the compactness of the SF re- 
gions. These range from very compact (<100 pc) nuclear star- 
bursts with almost no star-forming activity in the outer regions 
(type 1), to extended starbursts with relatively faint nuclei (type 
4), with type 2 and 3 as intermediate cases. In addition, they 
found a trend for galaxies with more compact SF region showing 
a higher star formation efficiency and hotter far-infrared color. 
They also found that the compactness of SF regions is weakly 
correlated with the galaxy morphology, with disturbed objects 
showing preferentially more concentrated SF. On the other hand, 
an appreciable fraction (~ 50%) of their galaxy sample was 
dominated by extended starbursts (type 4). The significant vari- 
ations in the degree of conc entration of the SB S F regions has 
been recently confirmed by McO uinn et al.l (120121) . 

In an attempt to characterize the SF spatial distribution for 
the SBs of our sample we derived estimates of thei r morpholog- 
ical cl ass and FIR colors. In particular, following Hattori et al. 
(12004 . SBs with strong tidal features and a single nucleus were 
classified as "mergers" (M), galaxy pairs with an overlapping 
disk or a connecting bridge were classified as "close pairs" (CP) 
if the projected separation is < 20 Kpc and galaxies that have a 
nearby (< 100 Kpc) companion at the same redshift were classi- 
fied as "pairs" (P). The remaining objects were classified as "sin- 
gle" (S). The classification of the SBs of our sample is listed in 
Tab.|6]along with the galaxy FIR colors, log/ 60 //ioo, log/2s//60- 

We attributed to each galaxy a compactness class on the ba- 
sis of its correlatio n with the FIR colors as shown in Fig. 4 of 
lHattori et al.l (|2004) that for the object of our sample corresponds 
to our Fig. [8] As it can be seen, we also confirmed their claim of 
a (weak) relation of FIR color and, as a consequence, compact- 
ness class with SB morphology. 

The next step is based on ISoifer etall ( 120001 12001 1) . For a 
number of SBs galaxies they plotted the MIR and NIR emis- 
sion curve of growth finding that in general the MIR emission 
is more concentrated, while only for few galaxies the MIR and 
NIR curves of growth show a similar trend. Actually we found 
that, to a first order approximation, the MIR emission profile of a 
given galaxy can be matched by NIR profile powered to an expo- 
nent a which ranges between 1 , when the two profiles are simi- 
lar, to 2, when the MIR emission is strongly concentrated. When 
we classify the sam e galaxies with the compactness criteria of 
lHattori et al.l (120041) . we found that (as expected) the galaxies 
with compact SF regions (type 1-2) are characterized by more 
concentrated MIR emission (a = 2.0 - 1.5, respectively), while 
galaxies with extended SF region (type 3-4) have similar MIR 
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Table 6: Morphological classification, FIR colors and com- 
pactness classificat ion for the SBs of our sample. Following 
lHattori et al. (2004), type 1 have SF region < 500 pc, type 2 < 1 
Kpc, type 3 > 1 Kpc and type 4 have extranuclear star formation. 



Galaxy 


Morp. 
class 


log/eo/Zioo 


log/25//eo 


Comp. 
class 


CGCGO 11-076 


S 


-0.20 


-0.89 


3 


CGCG043-099 


s 


-0.19 


-1.05 


3 


ESO148-IG002 


CP 


0.01 


-0.82 


2 


ESO239-IG002 


M 


-0.04 


-0.79 


2 


ESO244-G012 


CP 


-0.10 


-0.68 


2 


ESO264-G036 


s 


-0.34 


-0.96 


4 


ESO286-IG019 


M 


0.07 


-0.81 


2 


ESO440-IG058 


P 


-0.23 


-0.99 


3 


ESO507-G070 


s 


-0.08 


-1.21 


1 


IC1623A/B 


CP 


-0.14 


-0.80 


2 


IC2545 


M 


0.01 


-0.88 


2 


IC2810 


P 


-0.22 


-1.00 


3 


IC4687/6 


P 


-0.13 


-0.76 


2 


IRAS 12224-0624 


s 


-0.14 


-1.48 


1 


IRAS 14378-3651 


s 


-0.08 


-1.00 


1 


IRAS 16399-0937 


CP 


-0.24 


-0.87 


3 


IRAS 17207-0014 


M 


-0.05 


-1.36 


1 


IRAS 18090+0 130 


P 


-0.31 


-0.98 


4 


MCG-02-01-051/2 


P 


-0.14 


-0.79 


2 


MCG-03-04-014 


s 


-0.15 


-0.91 


3 


NGC0034 


M 


0.01 


-0.85 


2 


NGC0232 


P 


-0.23 


-0.89 


3 


NGC3110 


P 


-0.30 


-1.00 


4 


NGC5010 


s 


-0.33 


-0.85 


4 


NGC5331 


CP 


-0.29 


-1.00 


4 


NGC6240 


CP 


-0.06 


-0.81 


2 


NGC6926 


s 


-0.31 


-0.83 


4 


NGC7130 


p 


-0.19 


-0.89 


3 


NGC7592 


CP 


-0.12 


-0.92 


2 


NGC7674 


P 


-0.19 


-0.45 


2 
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Fig. 8: FIR color of the SBs of our sa mple. The differen t sym- 
bols identify the compactness class (Hatt ori et al.l 120041) while 
the label show the morphological type. 



vs. NIR profiles (a — 1.25 - 1.0, respectively). As a reference, 
we notice that in the typical case of NGC 6240, assuming a - I 
corresponds to locate 50% of the SFR within 1 .5 Kpc, whereas 
for a — 2 the same SFR fraction is enclosed within 500 pc. 

As a result of this discussion we have a prescription to es- 
timate the SF distribution based on the observed Lk map, and 
adopting a power index a appropriate for the compactness class 
of the given SB galaxy (Tab. [7]). 

3.1.5. Flow chart of the simulation 

Having defined all the ingredients of the simulation, we can now 
describe how this proceeds. The simulation flowchart can be 
summarized as follows: 

1 . for each galaxy of the sample, based on the estimated total 
SFR and adopted progenitor scenarios, we compute the ex- 
pected number of SNe per year; 

2. a time interval is chosen so that 100 SNe are expected to ex- 
plode in the given galaxy in that period. The time interval 
ends with the last observations of the galaxy. Given the ex- 
pected SN rates it is for all galaxies much longer than the 
duration of our monitoring campaign. The reason to simu- 
late 100 events is to avoid having to deal with fractional SN 
numbers for the different subtypes. We assign to each event 
a random epoch of explosion chosen within the defined time 
interval; 

3. each SN is assigned to a random explosion site inside the 
parent galaxy according to the SFR spatial distribution; 

4. a peak magnitude is also assigned to each SN, with a ran- 
dom value derived from the adopted SN luminosity function 
for the specific subtype. We also associate to the event an 
extinction value randomly extracted from a gaussian distri- 
bution whose mean value depends on the position of the SN, 
and cr - 1 /3 of the mean value; 

5. the apparent magnitude is thereafter calculated at the epochs 
of available observations, knowing the galaxy distance mod- 
ulus and SN epoch of explosion. This is compared with the 
search magnitude detection limit for the given position and, 
when brighter, the SN is added to the simulated discovery 
list. 

The process, iterated for all the galaxies of the sample, de- 
fines a single simulation run. Outcomes of the simulation are the 
expected number of SN discoveries, their types, magnitudes, ex- 
tinctions and positions inside the host galaxies. To explore the 
distribution of the outcomes from the random process, a com- 
plete experiment is made by collecting a minimum of a hundred 
single simulation runs. 

For the reference simulation we used theinput parameters 
summarized in Tab. [7] In the next section we will compare the 
prediction of the simulation with the current SN discoveries. 



4. Comparison between observed and expected SN 
discoveries 

As we outlined above, from a large number of MonteCarlo simu- 
lation runs we obtain the distribution of the expected SN discov- 
eries. This is shown in Fig. [9] where each bin of the histogram is 
the predicted probability of observing the specific number of SN 
discoveries whereas the dashed line marks the number of actual 
SNe discovered and the shaded area shows its 1-cr Poissonian 
uncertainty range. 
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Table 7: Input parameters for the reference simulation 



Ltir/S FR calibration 
IMF 
SFR distribution 

A v 

CC mass range 
la components 

SN luminosity function 



Kennicutt (1998) 
Salpeter 

for compactness class 1,2,3,4 
a = 2.0, 1.7, 1.25, 1.0 

oc S FR (max 30 mag) 

8 - 50M o 
30% prompty, 70% tardy 

LOSS (Li et al. 2011b) 



0.45 




16 17 18 

K magnitude 



0.20 



^0.15 



0.10 



IT} 
N 



2*0.05 
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Fig. 9: Histogram of the number of expected SNe from our 
MonteCarlo experiments. The dashed vertical line indicates the 
number of observed events and the grey area its 1 cr-Poissonian 
uncertainty. 



We found that with the adopted simulation scenario and input 
parameters we should have expected, on average, the discovery 
of 5.3 + 2.3 SNe. In 68% of the experiments (1-cr) the expected 
number is in the range 4-8 which is in excellent agreement with 
the observed number of 6 events. 

The prediction of the simulation is that almost all SNe are 
CC (5.1 SN CC vs. 0.2 SN la), though in 10% of the experiments 
at least one type la is found (that is what we have from the real 
SN search). 

The distributions of some of the expected and observed SN 
properties are compared in Fig. [TO] For the simulation, we show 
the distribution across a large number of experiments (line -only 
histogram) while the grey shaded histogram represents the actual 
observations. 

The top panel in Fig.[10]shows the distribution of the appar- 
ent magnitudes at the discovery. The good agreement between 
simulations and observations is a crucial consistency check of 
our estimates of the magnitude detection limit: if the discov- 
ered SNe were systematically fainter/brighter then expected, this 



(a) Expected (line-only) and observed (grey) K 
magnitude distribution at the discovery 




5 6 7 



A v mag 



(b) Expected (line-only) and observed (grey) ex- 
tinction distribution. In light grey we indicate the 
allowed range for the extinction of PSN2010 (see 
text). 



radial distance [Kpc] (NGC6240) 
10 5 3 1 0.50.3 0.1 




10 1 10 2 10 3 10 4 
K-band relative intensity 



(c) Radial distribution of expected (line-only) and 
detected (dark grey) SNe. For regular galaxies the 
surface brightness decrease monotonically with ra- 
dial distance (the upper axis shows this correspon- 
dence for one of the galaxy of our sample). In light 
grey we show the distribution of injected artificial 
SNe (see text). 
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would indicate, respectively, an underestimate/overestimate of 
the search detection efficiency. 

A comparison of the simulated vs observed extinction distri- 
bution is shown in the middle panel of Fig. [10] For the observed 
distribution the case of PSN2010 for which extinction is am- 
biguous is shown in light grey. Again the simulation is in good 
agreement with the observations. This argues in favor of the con- 
sistency of the input assumptions. The fact that, in our IR search 
we expect that most of detected SNe have low extinction (~ 75% 
with Ay < 1 mag) is a consequence of our assumption that the 
extinction is very high in the nucleus and rapidly decreases with 
the galaxy radius, following the same trend of the SFR. This 
does not means that extinguished SN are intrinsically rare but 
that they are confined to the galaxy nuclear regions where ex- 
tinction is extremely high even in the IR (see next). At the same 
time we can exclude the presence of a significant population of 
SNe with intermediate extinctions: we would easily detect them 
in our infrared search. 

Mattila& Meikiel feOOl) found an average value of Ay-30 mag 
for the extinction towards the SN remnants of M82. Confirming 
the presence of high extinction (about Ay ~ 15-45 mag) in the in- 
nermost 300 pc regions. On the other hand, Kanka re et al.1 (2008) 
found an host galaxy extinction of Ay ~16 mag for the SN 
2008cs, located at about 1.5 Kpc, relatively far from the galaxy 
nucleus. 

Finally, in the bottom panel of Fig. [10]we compare the distri- 
bution of locations inside the host galaxy for the expected (line- 
only) vs observed (dark grey) SNe. The different location are 
identified by the K band pixel counts: in general high counts 
occurs in the nuclear regions while low counts are in the out- 
skirts (we use pixel counts instead of radial distances because 
the latter is difficult to be defined for galaxies with irregular mor- 
phology or double nuclei. However a indicative correspondence 
from pixel count to radial distance is shown it the top axis of the 
figure for a galaxy with regular morphology. There is a mild in- 
dication of a deficiency of observed events in regions with high 
pixel counts. Taken to face value this may suggest a minor over- 
estimate of the detection magnitude limit in the nuclear regions. 
Given the poor statistics we cannot derive definite conclusions 
and therefore we will not elaborate further this issue. 

In the same figure we show also (in light grey) the distribu- 
tion of locations of the expected events for an ideal case where 
the magnitude detection limit in the nuclear regions is as deep 
as in the outskirts, and extinction is negligible. The experiment 
shows that the fraction of events that remains hidden to our 
search in the galaxy nuclear regions due to the combined effect 
of reduced search e fficiency and high extinction is very high, 
being about 60% (cf. lMattila et al.ll2012h . 

4. 1 . Uncertainties and effect of different assumptions 
Magnitude detection limit 

One of the main source of uncertainty for the simulation is re- 
lated to the estimate of the magnitude detection limit, magn m . For 
the reference simulation, we adopted as magum the mean value 
out of three artificial star experiments conducted for a number of 
selected positions inside the host galaxy (cf. l2.5b . The dispersion 
of measurements, that is the uncertainty on magn, n , is quite large 
with a typical range of ~ 0.5 mag but, in more extreme difficult 
cases, it can be as large as 2 mag. 

To test the propagation of this uncertainty, we performed 
MonteCarlo experiments assuming alternatively the lower and 
higher magum out of the three experiments. We found that the 



predicted number of SNe is respectively 6.2 + 2.5 and 4.7 + 2.2, 
that are +17% and -11% with respect to the numbers from the 
reference simulation. The fact that the error is significant is the 
reason why we spent a significant effort for a detailed estimate 
of the detection limit. 

SN Luminosity Function 

In the reference simulation we use a gaussian distribution for 
the SN luminosity f unction ( S N-LF) with a mea n value and dis- 
persio n taken from iLi et al ] (1201 lbl) . However, Hor iuchi et al.l 
based on a small sample of very nearby SNe, claimed 
that the faint end of the SN-LF is underestimated and SN CC 
fainter than mag -16 could made up to 50% of the dis- 
tributio n, to be compared wit h 20% of the sampl e of ILi et all 
d2011bl). On the othe r hand, iMattila et ail (1201 2|) argued that 
Horiuchi et al.l d201 ll) overestimated the fraction of intrisically 
faint CCSNe since they neglect the host galaxy extinction for 
their SN absolute magnitudes. We performed a MonteCarlo ex- 
periment adopting the Horiuchi's SN-LF and found that in this 
case the expected number of events would be low, only 3.3 on 
average. This is because most faint events are expected to fall 
below the search detection limit. The fact that the actual discov- 
eries are tw ice this number argues against a large fraction of faint 
SN-CC (cf. lBotticellaetal]l20Tl 

IMF and SN CC progenitor mass range 

The IMF enters both in the estimate of the number of SN pro- 
genitors and in the calibration of TIR luminosity in terms of SFR 
relation. However, the expected rate of CC SNe in our sample is 
virtually independent of the IMF slope. Indeed, for a given total 
mass of the parent stellar population, top heavy IMFs imply both 
a higher numb er of CC progenito rs as well as a larger luminos- 
ity. Following lDwek et al.l d201 ll) the number of CC progenitors 
per unit mass is kcc = 0.007, 0.01 1 and 0.039 M ~' respectively 
for a Salpeter, a Kroupa and a Starburst IMF , assuming that the 
progenitors range from 8 to 50 Mq. At the same time the total 
luminosity of a SB forming stars with a SFR of 1 Mq yr~' over a 
period of lOMyr (i.e. a 10 7 Mq stellar population) is 4.71 x 10 9 , 
7.33 x 10 9 and 2.55 x 10 10 L again for a Salpeter, a Kroupa and 
a Starburst IMF, respectively. The M/L ratio of such SB is then 
0.0021, 0.0014 and 0.0004 (solar units) for the three IMFs, and 
the expected number CC SNe originating from it is ^ 1 .5 every 
1O 5 L ' for all the three IMFs. Working out the numbers, it turns 
out that the SN CC rates from a population with a given Ltir is 
almost independent on the IMF, provided a consistent choice is 
made. 

Crucial is instead the assumption of the SN CC progenitor 
mass range, in particular the lower limit. Indeed if we adopted 
an upper limit of 100 Mq instead of the reference value of 50 
M Q the expected number of SNe would be 5.5 + 2.1, only ~5% 
higher then the reference simulation. On the other hand assum- 
ing a lower limit of IOMq (instead of 8 Mq ) results in an ex- 
pected number of SNe of 3.9+2. 1, which is ~30% lower than the 
expected rate obtained in the reference case. 

Extinction 

For the reference case we assumed that the extinction scales with 
the SFR with a maximum value corresponding to the SFR peak 
Ay = 30 mag. To test the uncertainty related to this assumption 
we made two different tests. In one experiment we maintained 
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the relation of Ay with SFR but taking, alternatively, a peak ex- 
tinction value A y = 10 and Ay = 100 mag. The experiment gave 
as expected number of SNe 5.5 + 2.7 and 4.7 + 2.3, respectively. 
In the second experiment we assume that the extinction is con- 
stant through the galaxy and is Ay = 3.0 mag. In this case the 
expected number is 5.3 + 2.3 identical to the value of the refer- 
ence simulation. 

The conclusion is that the uncertainty on the extinction does 
not affect significantly the simulation or, conversely, that our ex- 
periment we cannot probe the extinction distribution. 

Star Formation Distribution 

The spatial distribution of SFR is an important, and the most 
uncertain, ingredient of the simulation. For instance, if we as- 
sume that the SFR is confined in the very inner regions, say in 
the inner 3 - 500 pc, the resulting SNe will remain unaccessi- 
ble to our search. On the other hand, the fact that in some SBs 
the SFR is extended h as been confirmed by different studies (eg. 
McO uinn et al ] l20l2h . not to mention that many of the SNe we 
have discovered are at significant radial distances (cf. Tab. [3}. 

As we described in Sect. 13.1.41 as proxy of the SFR distri- 
bution we use L a K where a range from 1 to 2 depending on the 
galaxy compactness class (Tab. To test for the uncertainties 
of this assumption we performed two simulations assuming that 
for all galaxy a is either 1 or 2. We obtained in the first case 
an expected rate of 8.8 + 3.0 and in the second case a value of 
3.0 + 1.7. The latter occurs because when the SFR is more con- 
centrated, a large number of SNe remain hidden to our search 
due to the low search detection efficiency in the nuclear regions. 

The conclusion is that the uncertainty in the adopted SF dis- 
tribution propagates with an error of ~ 50% on the expected SN 
number. We may consider that the actual good match of observa- 
tions with the reference simulation argues in favor of the adopted 
prescription. 

5. Summary and Conclusions 

We have presented the analysis of an infrared SN search in a 
sample of 30 nearby SB galaxies, conducted between 2009 and 
2011, with the goal to verify the link between star formation 
and SN rate. During our search we collected in total about 240 
observations discovering 6 SNe, 4 of them with spectroscopic 
confirmation. 

How does this number compares with the expectation ? 

Answering this question requires a detailed characterization 
of the SN search detection efficiency, the galaxy properties (in 
particular SF rate and spatial distribution) and the SN properties 
and progenitor scenarios. We included all these ingredients in 
a MonteCarlo simulation tool that, allowing for the stochastic 
nature of SN events, can be used to explore the distribution of 
the expected SN number and properties. 

First of all, we may remark that by itself the number of de- 
tected SNe is a proof of the high SFR in SBs. In fact if we com- 
pute the expected number of SNe in our surve y based on th e 
average SN rate per unit B luminosity or mass dLi et al.ll2.Ql lab , 
we would predict the discovery of 0.5 events (or more precisely, 
50% of the simulation predict the discovery of one event and 
none is expected in the other 50%). The observed number is one 
order of magnitude higher, which is consistent with the fact that 
the TIR emission of SBs is about ten times higher than for nor- 
mal SF galaxies with the same B luminosity. Indeed, it is well- 
known that the TIR luminosity is an excellent tracer for SFR, in 
particular in SBs. 



When we adopt the SFR from Ljir as input for the 
MonteCarlo experiment, we find that the expected number of 
SNe in our search is 5.3 + 2.3, SNe in excellent agreement with 
observations. In most cases we predict that only SN CC should 
be discovered while in the actual search we did detect one type la 
SN. Given that there is a sizable fraction of experiments (10%) 
when this is predicted to occur we do not elaborate further this 
issue. Also, allowing for the low statistics, we find an excellent 
agreement between the predicted and observed SN properties, 
namely apparent magnitude at discovery, extinction and location 
inside the host galaxies. 

We performed a number of tests to verify the dependence 
of the simulation outcomes from the input parameters. For the 
SN search characterization we show that an accurate estimate of 
the magnitude limit for SN detection is crucial. This is why we 
spend a considerable effort in artificial star experiments (possi- 
bly the single most expensive task of our project). For the galaxy 
characterization the most uncertain input is the SF spatial dis- 
tribution. With some creativity, we devised a prescription that 
seems to work, but it is certain that this is a place for improve- 
ments when new, high resolution SB maps will become avail- 
able. Instead, we found that our results are not sensitive to the 
uncertainty on the amount of extinction because where extinc- 
tion is very high (the dense SB regions) our search is limited by 
the bright magnitude detection limit. SNe in these regions re- 
main hidden to our search almost independently on the amount 
of extinction. Based on our simulation we estimated that the frac- 
tion of hidden SNe is very significant, that is ~ 60% with an up- 
per limit of 75% if we account for the poissonian uncertainties 
in the number of detected events. Finally, for the SN progenitor 
scenarios the larger uncertainty is the lower limit of the progen- 
itor mass range. If we adopt a lower limit M^ c = 10 Mq instead 
of 8 M0 as in the reference simulation, the expected number of 
SNe would be 30% lower than observed. 

Our results appea r in good agreement with those of previ- 
ous similar searches dMannucci et al.l 120031; ICresci et ail 120071: 
Mani la et al.ll2007all2012l cf. Sect. PP. In broad terms, the over- 
all conclusion of all these studies can be expressed as follows: 
the number of (CC) SNe found in SBs galaxies is consistent with 
that predicted from the high SFR (and the canonical mass range 
for the progenitors) when we recognize that a major fraction 
of the events r emains hidde n in th e unaccessible SB regions. 
As stressed by iMattila et al.1 (120121) . this has important conse- 
quences for the use of SN CC as probe of the cosmic SFR, be- 
cause Jfie_jiracJj£nof_SBs^ with redshifts 
(cf . lMelinderetaLll2012l: iDahlen et al.ll2012l) 

While continuing to search for SNe in SBs, in optical and 
infrared, can certainly help to improve the still low statistics, one 
may argue at this point for a change of strategy. 

In this respect good example is the attempt to reveal some 
of the hidden SN CC through infrared SN searches which ex- 
ploits adaptive optics at large telescopes, eg. Gemini or VLT. 
The results are encouraging with the discovery of two SNe with 
very high ex tinction, namely SN 2004ip with Ay between 5 
and 40 mag (IMattila et alJl2007bh and SN 2008cs with Ay 16 
mag (IKankare et alJl2008l) . though we may notice that both ob- 
jects were too faint for spectroscopic confirmation. Other two 
SNe were discovered very close to the galaxy nucleus, namely 
SN 201 0cu at a radial distan ce of 180 pc and SN 2011 hi at 
380 pc (IKankare et al.l 120121) . though in these cases the low 
extinction suggests that the low radial distance is a projec- 
tion effect (also in these cases no spectroscopic classification 
wa s obtained). The extinction tow ards SN 201 lhi was revised 
by iRomero-Cafiizales et al.l d2012l) using Gemini-N data. They 
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demonstrate that this is most likely a SN IIP with A v of 5-7 mag. 
Because of the need to monitor one galaxy at the time and to ac- 
cess heavily subscribed large telescopes, this approach will not 
result in large statistics though even a few events may be molstly 
valuable to explore the very obscured nuclear regions. 

On the other hand, a new opportunity that should be explored 
is the piggy-back on wide field extragalactic surveys of the next 
generation infrared facilities, in particular EUCLID. This would 
allow for the first time to perform IR SN searches on large sam- 
ple of galaxies exploring a range of SF activity and, by monitor- 
ing galaxies at different redshifts, probe the cosmic evolution. 
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